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Abstract 

The neutron star crust is rigid material floating on a neutron-proton 
liquid core. As the star's spin rate slows, the changing stellar shape 
stresses the crust and causes fractures. These starquakes may trigger 
pulsar glitches as well as the jumps in spin-down rate that are observed 
to persist after some glitches. Earlier studies found that starquakes in 
spinning-down neutron stars push matter toward the magnetic poles, 
causing temporary misalignment of the star's spin and angular momen- 
tum. After the star relaxes to a new equilibrium orientation, the mag- 
netic poles are closer to the equator, and the magnetic braking torque is 
increased. The magnitude and sign of the predicted torque changes are 
in agreement with the observed persistent spin-down offsets. Here we 
examine the relaxation processes by which the new equilibrium orienta- 
tion is reached. We find that the neutron superfiuid in the inner crust 
slows as the star's spin realigns with the angular momentum, causing 
the crust to spin more rapidly. For plausible parameters the time scale 
and the magnitude of the crust's spin up agree with the giant glitches 
in the Vela and other pulsars. 



INTRODUCTION 

Stresses in the crust of a neutron star could produce starquakes that 
affect the star's spin evolution and generate high-energy emission. As 
the star's spin rate increases or decreases, changes in the equilibrium 
shape of the star and the differential rotation between the crust and the 
interior neutron superfiuid generate stress [1]. In "magnetars", decay 
of the superstrong field {B ^ 10^^ G) could break the crust and drive 
episodes of intense gamma-ray emission [2]. Recent studies showed that 
starquakes can change the magnetic spin-down torque acting on the 
star [3, 4]. Starquakes in slowing neutron stars drive matter toward 
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the magnetic poles, distort the star's shape, and excite precession. As 
the precession damps, the star relaxes to a new rotational state with 
the magnetic poles closer to the equator. The new magnetic orientation 
enhances the braking torque on the star and may provide an explanation 
for the observed increases in the spin-down torque following glitches in 
the Crab pulsar, PSR1830-08 and PSR0355+54. 

Here we investigate the physical processes that allow the star to relax 
to its post-star quake equilibrium. The most important processes are 
the coupling between the liquid core and solid crust and the creep of 
neutron superfluid vortices in the inner crust of the star. We find that the 
changes produced by large starquakes can trigger catastrophic unpinning 
of neutron superfluid vortex lines in the star's inner crust [5]. As vortices 
move, the inner crust superfluid rapidly settles to a state of lower angular 
momentum, while exerting a spin-up torque on the crust. Our estimates 
show that starquake-triggered events may explain giant pulsar glitches 
as well as the persistent spin-down offsets. The next section summarizes 
earlier work on starquakes and spin-down off'sets, and the following one 
describes the post-starquake spin relaxation and glitches. 

STARQUAKES AND SPIN-DOWN OFFSETS 

The crust of a spinning neutron star is oblate with an equatorial bulge. 
The moment of inertia of the bulge is leb ~ laustR^^'^ / C^GM) ~ 4 x 
10"^/crust^^2 where R ~ lO'^ cm is the stellar radius, O = IOOO2 rad s is 
the star's spin frequency, and Icrust is the characteristic moment of inertia 
of the crust. The crust contains about 1% of the star's total moment 
of inertia /total- As the star spins down, the equatorial circumference 
shrinks and the polar radius grows. Because the crust is solid, strain 
develops as the star's shape changes. As sketched in Figure 1, the strain 
in the crust is relieved along starquake faults that form at an angle to the 
star's equator [3, 4]. Matter slides along these faults to higher latitudes, 
and magnetic stresses favor those faults that direct matter along field 
lines toward the magnetic poles. 

An important result of the earlier studies is that starquakes shift 
the stellar matter asymmetrically, creating excess moment of inertia 61 
about an axis different from any of the pre-starquake principal axes. This 
distortion changes the orientation of the principal inertial axis by an an- 
gle Aa ~ 6I/Ieb ~ 2.5 X lO'^n^^ 61 /Icrust = 2.5 X 10-^^2^6-7, where 
61 /Icrust = 10"'' (5-7 [3]. The distortion parameter 6-7 characterizes the 
size of the starquake. 

When the principal axis of inertia of the crust is not aligned with the 
star's angular momentum, the star precesses and wobbles. Eventually 
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the star relaxes to a new equilibrium in which the axis and angular 
momentum are again realigned, and the magnetic pole is shifted by 
the angle Aa closer to the rotational equator. In some mechanisms 
for pulsar spin down, such as the magnetic dipole braking model, this 
angular shift increases the torque on the star, producing a long-lasting 
increase in the spin-down rate. The persistent spin-down offsets observed 
in the Crab pulsar can be explained by this mechanism if Aa ~ 10~^ [6] 
corresponding to 5-7 ~ 1.6 ( 1^2 — 2 for the Crab pulsar). 




Aa 



Figure 1 Starquakes re- 
lieve the stress that builds 
as the a neutron star's spin 
slows. Matter can slide 
to higher latitudes along 
faults F or F' . Mag- 
netic stresses favor faults 
such as F that move mat- 
ter toward the magnetic 
poles. The accumulated 
matter, shown as snow- 
capped peaks, shift the 
principal axis of inertia by 
an angle Aa relative to the 
star. 

THE ORIGIN OF GLITCHES 

Starquake-induced asymmetry in the stellar crust excites precession. 
If the neutron star crust behaved as an isolated rigid body, it would pre- 
cess or wobble at a frequency Qyj ~ {leb/ Icrust)^ and the angle between 
the angular velocity and the angular momentum would be ^ (Aa)^. The 
spin behavior of a realistic neutron star is more complicated than this 
for several reasons. First, the pinning of superfluid vortices in the inner 
crust acts to stabilize the spin of the crust. Second, the crust and core 
of the star are not strongly coupled on the precession time scale. Third, 
the tilting of the crust accelerates vortex creep. 

Pinned superfluid. In the inner crust of the star the neutron su- 
perfluid vortex lines may pin to the nuclei in the solid crust. The rota- 
tion of the superfluid is determined by the location of the vortex lines. 
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and, as long as the vortex lines remained pinned, the superfluid veloc- 
ity field cannot change. The gyroscopic action of the pinned super- 
fluid works with the equatorial bulge to further stabilize the star [7]. 
The moment of inertia /pinned of the pinned superfluid is comparable 
to that of the crust and much larger than the moment of inertia of 
the equatorial bulge; /pinned ^ 0.01/totai ^ leb- The pinned superfluid 
decreases the equilibrium tilt of the star by a factor /e6//pinncd to an 
angle Aap ~ (^///pinned 10~^5_7 -C Aa. This is the tilt angle of the 
star with completely pinned superfluid after the precession has damped 
(the tilt immediately after the starquakc is 

The precession frequency is proportional to /pinned inversely pro- 
portional to moment of inertia coupled to the crust. If the spin of the 
core and the crust are tightly coupled, the star's precession frequency is 

= (/pinned //total — ^^2 I'ad S ^. 

Coupling between the crust and the core. Changes in the crust's 
motion are communicated to the core by MHD-like waves. If the protons 

in the neutron star core form a type II superconductor, as expected, the 
magnetic field is confined to thin tubes of flux irhc/e ~ 2 x 10~^ G 
cm~^, with characteristic dimensions A ~ 50 fm and field strengths 
~ 10^^ G. Signals travel from the crust to the core on an Alfven 

time tcoupic ~ [AiTp)^/'^ R/ {B B^)^/"^ ~ 4:B^2^'^ s, for an average field of 
B = 10^^i?i2 G and a density of p ~ 10^^ g cm~^ [8]. Because of the 
imperfect coupling between the crust and the core, the precession is 
damped on a time scale [9] 

idamp — T^icouple ~ J'"^'^' ^couple ~ 4005^2^''^ S. (1) 
^'lo -'pinned 

The energy of the regular motion of precession is converted into irreg- 
ular fiuid motions in the core. The irregular motions of the superfiuid 
neutrons and superconducting protons are then converted into thermal 
energy by processes such as electron scattering from vortex lines [10] or 
flux tube- vortex line interactions [11]. Since the core and crust arc not 
strongly coupled on the precession time scale, the effective moment of 
inertia of the precessing material is less, and the precession frequency 
may be higher than we have used. The damping time scale would be 
correspondingly reduced. 

After the precession has damped, the star's crust has tilted by Aap. 
As the pinned vortex lines move to align with the crust's rotation vector, 
the tilt angle grows until it equals Aa. We now examine the processes by 
which the vortex lines move and the effects of this motion on the star's 
spin behavior. 
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Figure 2 The shift of the 
angular velocity from Q, to 
D! increases the velocity 
lag between the superfluid 
and the solid crust in the 
upper right and lower left 
sections of the star's inner 
crust in this sketch. In 
these sections the vortex 
creep accelerates and the 
vortex lines migrate out- 
ward and align with the 
new direction of the angu- 
lar velocity. 

Motion of vortex lines. In a steady state, the angular velocity of 

the superfluid closely follows that of the solid crust. As the superfluid 
rotation slows, the vortex lines move radially with a steady velocity 
Vss — —R(i/{2Q) = i?/(4tage), where the spin-down age of the pulsar 
is tage = ^/i^tl) [14]. For the Crab pulsar Vss — 5 x 10~^ cm s"-*^ and 
for the Vela pulsar v^s — 7 x 10"''' cm s~^. If the vortex lines move at 
Vss following a starquake, they would align with the new direction of 
the crust's principal axis in a time taiign ^ ^dR/vss ~ lO^^^-yfi^^tage, 
years for the Crab pulsar and decades for the Vela pulsar. The vortex 
velocity can be greatly enhanced by the perturbations produced by the 
tilt of the spin axis.^ After the star's spin axis tilts relative to the pre- 
starquake orientation, some parts of the star are further from the new 
spin axis; this shift can be as large as AapR. For these regions, the 
linear velocity from the star's spin increases by AapR^l. The superfluid 
velocity, on the other hand, changes much less. If the neutron vortex 
lines in the star vortex were completely pinned, the superfluid velocity in 
the star's frame would remain unchanged. Even though the superfluid 
in the core is not pinned, the vortex lines cannot move through the 
crust, so the crust superfluid velocity remains fixed. The velocity lag 
vs between the superfluid and the solid crust thus changes by as much 
as Avs ~ AapR^} ~ 10^6^7 cms~^ in parts of the star; see Figure 2. 
Velocity differences of this magnitude may have dramatic effects on the 




'^Heat produced during a starqucike can also increase the vortex velocity [12]. 
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vortex velocity, large enough to produce the crustal spin-ups associated 
with glitches. 

Pinned vortex lines in the inner crust can creep outward through 
thermal activation of vortex segments from their pinning barrier; a line 
segment unpins from one configuration and migrates to new sites where 
it repins [13, 14]. The general form of the creep velocity is Vcreep oc 
exp[—A/kT], where A is the minimum activation energy for a segment of 
the vortex line to unpin [15]. The activation energy depends sensitively 
on the lag velocity vs- If the vortex lines are moving at their steady state 
rate Vss and the lag velocity suddenly jumps by Avs, the new creep rate 
is 

^^creep-^^..exp^ dt-^ kT \ ' ^ ^ 

The activation energy of a vortex line is characterized by the maximum 
pinning force Fpin between a nucleus and a vortex line, the range of 
the pinning force rpin and the distance between pinning sites ^pin- The 
maximum lag velocity between the superfluid and the crust that the 
pinning forces can support is i^max = Fpin/ Ps ^'^pin, where k = nh/ (m„) = 
2.0 X 10~^ cm^ s~^ and ps is the mass density of the neutron superfluid; 
Ps ~ 10^^ g cm~^ is the characteristic density of much of the inner 
crust superfluid. The activation energy required for a vortex line to 
unpin depends on the stiffness of the vortex line; if the line is flexible, 
it can unpin from one nucleus at a time, whereas if it is stiff, the line 
must unpin from many nuclei simultaneously [15]. The parameter r ~ 
0.4psK^rpin/(Fpi„£p;ii) characterizes the vortex line's stiffness [15]; for 
the conditions in much of the crusts of the Crab and Vela pulsars r > 1 
and vortex lines move by unpinning from many sites simultaneously. In 
this limit the appropriate expression for the activation energy is ^ ~ 
6.8FpinrpinT^/^(l — vs/vj^ax)^^'^, and it's derivative is 

^ ^ _4.8^PHi!>^ . -3.0 ( ^1 . (3) 

dvs 

In obtaining this equation we set (1 — vs/vraax)^^^ ~ 0.6, which is a 
characteristic value for a variety of stellar models [14]. 

We ca,n use the observRitions of the Vela pulsar to estimate I'max* 

This 

pulsar exhibited a string of 12 nearly evenly spaced glitches separated 
with an average interval of tint = 2.3 years [16]. The regularity of these 
glitches suggests that the inner criist superfluid remains pinned between 
glitches until the velocity lag between the solid crust and the superfluid 
approaches Umax- With this interpretation, a lower limit to the critical 

lag velocity is Vmax~ I^^VelaKint-R = 7.0 X 10^ cm S~^. The value of Vmax 
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exceeds this limit because the inner crust superfiuid might not relax to 
zero lag velocity at each glitch. For our estimates we use t^max = lO^vs cm 
s~^. For ~ 1 the pinning force is -FpinPs'^^pin^max — 13 (£pin/100fm) 
keV fm^^ at a superfiuid density of ps = 10^^ g cm~^. This pinning force 
is much smaller than that obtained by recent microscopic calculations 
of vortex-nuclear interactions [17] assuming a perfect crystal, but it is 
larger than the average pining force estimated for an amorphous crust 
[18]. 

Taking pg = 10^^ g cm ^ and interaction distance rpin = 10 fm [19], 
the change in the activation energy in parts of the crust after a starquake 
is 

^Avs ~ -0.053t;r^/^ ( -^^] keV ~ -52 6-7^2v7^^'^ keV. (4) 
dvs V cms~^/ 

The time scale for the post-starquake relaxation of the vortex lines is 
^reiax ~ 'ActR/vcTeep'i for kT in keV, this gives 



^relax ^ tSt SXp 



52 5-7n2 



(5) 



The vortex relaxation time is very sensitive to the magnitude of the 
quake- induced shape change. For example, taking fs = 1, we find that 
for the Crab pulsar {kT ~ 20 keV, tage ~ 10^ yr, Bu ~ 4) the relax- 
ation time is less than the damping time if the distortion parameter is 
(^-7~8.8, but it is more than lOOOtdamp if (5-7^4.3. The corresponding 
values for the distortion parameter for the Vela pulsar (kT ~ 8 keV, 
^agc 10^ yr, Bi2 ~ 8.8) are and S-7'^1.2. For each pulsar, 

there is a critical value for the distortion parameter S-7. Star quakes 
that produce distortions above this threshold trigger rapid vortex mo- 
tions while smaller events generate only gradual changes. 

Spin jumps. In the large events, the rapid outward motion of the 
vortex lines produces a slowing of the superfiuid and a corresponding 
spin up of the crust. In the regions of the inner crust where Avs is posi- 
tive the vortex lines rapidly creep a distance AaR 

The superfiuid in the affected regions of the inner crust will now spin 
more slowly by corresponding amount: Ai^s/^ ^ 2.5 x 10^'^$7^^(5_7. 
The angular momentum lost by superfiuid is imparted to the crust and 
the core which is strongly coupled to it [10], giving /total ^^crust = 
— /mipijincd'^^^s) where /unpinned is the moment of inertia of the region 
of rapid vortex creep. Since only 1/2 of the crust has Avs > 0, we have 

-^unpinned ^ 0.5/crust &ud 

^^crust ^ _ -^crust A^lg ^ ^ 3 x 10~^*^~^ (6) 
^ 2/total ^ ^2 
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The rapid creep of vortex lines following a large starquake could ex- 
plain the giant glitches with AJlcrust/^ ~ 2 x 10""^ observed in the Vela 
and other pulsars. 

In the above estimate we assume that the lag between the superfluid 
and the crust is large enough to supply the needed angular momentum 
to the crust. As the vortex lines creep through a distance ^ AaR the 
local superfluid slows by ~ AaRO,. If the pre-starquake lag velocity 
vg were less than this value, the vortex creep would stop before the 
vortex lines move the full distance ~ AaR. A necessary condition for 
the above estimate to be valid is that vs ^ AaQR. Since i^max > ^<5 

we have v^^2.55-'i'Q2^ ■ '^^^ limiting factor in the size of a starquake- 
induced pulsar glitch may be the pre-glitch lag velocity. The reason for 
the large difference between the size of the Crab glitches and the giant 
glitches of the Vela pulsar may be that in the Crab pulsar, with its higher 
internal temperature, vortex creep between glitches limits the build up 
of a sufficiently large vs- 

SUMMARY 

Starquakes tilt the principal axis of inertia of a neutron star. As 
the star relaxes to its new equilibrium orientation, neutron superfluid 
vortex lines migrate outward, spinning up the rest of the star. The 
spin ups from large quakes may explain the giant glitches observed in 
isolated pulsars. The change in the direction of the magnetic axis may 
increase the spin-down torque, as observed in the Crab and other pulsars. 
Detailed calculations of the post-starquake relaxation for both large and 
smaller events may yield distinctive timing signatures to compare with 
observations. 
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